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Suitable shaping (in particular, flattening and broadening) of the laser beam has recently been 
proposed as an effective device to reduce internal (mirror) thermal noise in advanced gravitational 
wave interferometric detectors. Based on some recently published analytic approximations (valid 
in the infinite-test-mass limit) for the Brownian and thermoelastic mirror noises in the presence of 
arbitrary-shaped beams, this paper addresses certain preliminary issues related to the optimal beam- 
shaping problem. In particular, with specific reference to the Laser Interferometer Gravitational- 
wave Observatory (LIGO) experiment, absolute and realistic lower-bounds for the various thermal 
noise constituents are obtained and compared with the current status (Gaussian beams) and trends 
("mesa" beams), indicating fairly ample margins for further reduction. In this framework, the ef- 
fective dimension of the related optimization problem, and its relationship to the critical design pa- 
rameters are identified, physical-feasibility and model-consistency issues are considered, and possible 
additional requirements and/or prior information exploitable to drive the subsequent optimization 
process are highlighted. 

PACS numbers: 04.80.Nn, 07.60.Ly, 41.85.Ct, 42.55.-f 
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I. INTRODUCTION 



In all currently operating (and possibly future) inter- 
ferometric gravitational wave detectors, the overall limit 
sensitivity of the instrument is bounded by the noise 
floor, which, in the most interesting observational fre- 
quency band (30-300 Hz), is dominated by thermal noises 
in the substrate and in the high-reflectivity coating of 
the test masses. With particular reference to the Laser 
Interferometer Gravitational- wave Observatory (LIGO) 
experiment [l[ , an introductory discussion of the various 
noise components can be found in (2j, and a numerical 
code for computing the noise budget is available from 
Q . Toward the development of second-generation detec- 
tors, such as Adv-LIGO [U, the quest for increasing the 
event rate in the observational band has motivated the 
exploration of various techniques for reducing the mir- 
ror thermal noise. With specific reference to the coating 
Brownian noise (dominant in the current baseline design 
featuring fused-silica test-masses) , use of improved (low- 
mechanical- loss) materials Q , geometric optimization of 
the coating design [f|, and flat-top (commonly referred 
to as "mesa") beams @, § seem the most promising. 
The latter option, intuitively motivated by the poten- 
tial capability of a mesa beam (MB) of better averag- 
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ing the thermally-induced mirror surface fluctuations as 
compared to a standard Gaussian beam (GB), has been 
numerically proved to yield significant reductions in the 
overall thermal noise [e], [l(| , and has led to the develop- 
ment of a cavity prototype with non-spherical "Mexican 
hat" (MH) profile mirrors [TJ E2] • Alternative (nearly- 
concentric [13|], nearly-spheroidal [3, EH, EH) designs 
have been subsequently proposed to cope with the in- 
herent tilt-instability of the originally-conceived nearly- 
flat configuration. Also, use of higher-order modes in 
standard spherical cavities has been shown to provide, 
in principle, comparable reductions without the need of 
changing the mirror profile [ItJ , but its practical feasibil- 
ity still remains to be assessed. 

The method utilized in [lCj to compute the coat- 
ing and substrate thermal noises relies on a finite-test- 
mass (FTM) computationally-intensive numerical anal- 
ysis based on the approach in [IH, EH- More recently 
pfl [2l| . a general though simple formula has been de- 
rived in the infinite-test-mass (ITM) limit, which allows 
the computation of the above noises for arbitrary-shaped 
beams. This approximation has been validated and cali- 
brated in [2l[ against the FTM numerical solutions (see 
also the discussion in Sec. IIII B 31) . In view of its remark- 
ably simple form, in terms of spectral integral function- 
als of the beam intensity profile, it appears suggestive to 
exploit it for addressing the optimal beam-shaping prob- 
lem, i.e., finding the beam profile that minimizes a given 
thermal noise constituent. In a step- by-step approach, 
acknowledging the formal and computational complexity 
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of the arising optimization problem, this paper addresses 
some key preliminary issues. In particular, emphasis is 
placed on the a priori deduction of absolute and realistic 
lower-bounds for the various thermal noise constituents, 
the identification of the effective dimension of the prob- 
lem, and how this depends on the critical design parame- 
ters, and the gathering of additional requirements and/or 
prior information to be utilized in the actual optimization 
problem. 

Accordingly, this paper is organized as follows. In Sect. 
HI1 the problem geometry, formulation, notation and 
strategy are outlined, with a compact review of the rele- 
vant background theory (ITM approximation). In Sect. 
Mil under the idealized assumption of zero diffraction-loss 
{compact spatial support) beam profiles, absolute lower 
bounds for the noise constituents, as well as the corre- 
sponding beam profiles over the mirror, are obtained in 
analytic form, by solving a straightforward variational 
problem. Subsequently, a key physical-feasibility con- 
straint (related to the finite spatial bandwidth of the cav- 
ity eigenmodes) is taken into account by approximating 
the above compact-support optimal profiles in a suitable 
L 2 functional subspace, whose dimension is fixed by the 
diffraction-loss constraint. This results in more realis- 
tic tighter bounds. In this framework, the role of the 
number of electromagnetic degrees of freedom [22j of the 
cavity in setting the effective dimension of the optimiza- 
tion problem is highlighted. In Sect. IIV1 the obtained 
absolute and realistic lower bounds for the considered 
noise components are compared to the levels currently 
achievable using GB and MB profiles. Moreover, some 
model-consistency issues are discussed in order to assess 
the practical relevance of the results. Finally, in Sect. [Vj 
conclusions and recommendations are provided. 



II. PROBLEM STATEMENT 



A. Geometry 



Referring to the problem geometry illustrated in Fig. 
[TJ we consider a standard Fabry-Perot optical cavity with 
two identical, symmetric, nearly-flat (nonspherical) mir- 
rors of radius a laid on cylindrical test masses, separated 
by a distance L (see Fig. [TJi). The mirror (axisymmetric) 
departure from flatness is described by h(r), with r de- 
noting the radial coordinate in the mirror plane (see Fig. 
[T}d). In what follows, attention is focused on the axisym- 
metric (i.e., ^-independent) eigenmode field distribution 
$(r) on the mirror, with implicit exp(iwi) time-harmonic 
dependence. Note that, in view of the duality relations 
expounded in [2j|, the results derived hereafter apply to 
the nearly- concentric case too 



(a) 

t--h(r) i-h(r 
nearly-flat mirrors 



L 

test masses 




FIG. 1: Problem schematic: A perfectly symmetric Fabry- 
Perot optical cavity composed of two nearly-flat mirrors (with 
profile h(r)) attached on cylindrical test masses of radius a 
separated by a distance L along the z-axis. (a): Side view, 
(b): Front view. 



B. Background: Infinite- Test-Mass 
Approximations 

In the ITM approximation [2(| Hl[ , and in the low fre- 
quency limit of interest for gravitational-wave interfer- 
ometers, the power spectral densities of the main coating 
and substrate thermal noise constituents of interest can 
be written as 



S = C 



poo 2 
J K q+1 |W [|$| 2 J («)} dK 



(1) 



where C is a noise-type- and frequency-dependent factor 
(irrelevant for all further developments), q is a noise- type- 
dependent scaling exponent (see Table [J), <I>(r) is the 
axisymmetric eigenmode field distribution on the mirror, 
and 



(2) 



denotes the Hankel-transform (HT) operator. Here and 
henceforth, J m denotes an mth-order Bessel function of 
the first kind 0, Sec. 9.1]. The (axisymmetric) field 
distribution <f>(r) satisfies the eigenvalue equation [25| 



7<E>(r) 



K(r, r')$(r')r'dr', 



(3) 



where 7 denotes the half-roundtrip eigenvalue, and the 
kernel is given by 



x exp < ik 



h{r) + h{r') 



(r 2 +r' 2 ) 



2L 



■(4) 



with k — 2n/X denoting the free-space wavenumber (A 
being the wavelength). Equations ([3]) and (j4]) can be 
recognized as a mapping between a mirror profile h(r) 
and a set 0(/i) = {[j m , $ m ], m = 1, 2, . . . } of eigenstates. 
Here and henceforth, unless otherwise specified, the field 
distribution on the mirror is assumed to be normalized 
as follows 



|$(r 



rdr = 1. 



(5) 
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TABLE I: Thermal noise constituents of interest and corre- 
sponding scaling exponents (cf. JTJ). 
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with Ct denoting a design limiting value. For Adv- 
LIGO, the reference figure is Ct — lppm (10~ 6 ). The 
diffraction-loss constraint singles out a subset flc(h) C 
Q(h) of admissible eigenmodes. 



C. Formulation and Notation 

It is expedient to recast the problem into a canonical 
form by introducing the scaled variables 



K 



an, 



and the scaled field distribution 

4>{f) = a$(fa). 



(7) 



(8) 



Here and henceforth, the overbar denotes scaled quan- 
tities. The noise functional in |T|) can accordingly be 
rewritten as 



S = 



C 



■S 



(9) 



where 
S 



R q+1 \n [|</>| 2 ] (s)} 2 ds, (io) 



thereby explicitly factoring out the a~^ 9+2 - ) scaling law 
predicted by the ITM approximation [2(1 HH • In what 
follows, we focus on the scaled noise functional in ([TO]) , 
which essentially accounts for the beam-shaping effects. 
Unless strictly needed, the explicit dependence on \<j)\ 2 
and q will be omitted for simplicity of notation. The 
scaled field distribution <j)(f) in §8§ satisfies the scaled 
version of the eigenproblcm in ([3]), which can be conve- 
niently recast as 

70(f) = iitNd exp [— iV(f)} Hi [exp (— iV) 4>] (TriYof), 

(11) 

where 7 = 7 exp (ikL), 



ni[m)= I 

Jo 



F(oj mcdc 



(12) 



denotes the [0, 1] interval-windowed HT operator, and 



V(f) = kh{af) 



TlN D f 2 



is a mirror-profile-dependent phase function, with 



Nr 



2N F = 



2a^ 
AL 



(13) 



(14) 



denoting twice the Fresnel number Np of the optical cav- 
ity [HI . Following , we shall refer to N D in (HU) as the 
number of electromagnetic degrees of freedom [39(, whose 
relevance will be illustrated later on (see Sec. IIII Bj) . In 
the following we shall always assume the eigenfunctions 
as normalized, viz., 



|(/>(r)| 2 rdr 



= 1, 



(15) 



with ||-|| denoting the usual VL , (cylindrical) Hilbcrt 
norm. Accordingly, we shall write the diffraction-loss 
constraint as 



C\6] = 



\(f){f)\ 2 fdf = 1 - |7| 2 < C T - (16) 



D. The Optimization Problem 

The optimization problem of interest consists of mini- 
mizing the scaled noise functional in (|10[) . acting on the 
mirror profile h(f), i.e., in finding the special mirror pro- 
file h(f) (within a suitable functional class, e.g., C°°) for 
which 
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(17) 



takes on its smallest value, Q,c(h) denoting the subset of 
eigenmodes obeying the diffraction-loss constraint (|16[) . 
The minimization of (fT7| . subject to p^|) . represents a 
formidable optimization problem, whose well-posedness 
(i.e., existence and uniqueness of the solution, and its 
continuous dependence on data) cannot be taken for 
granted, with the consequent ill- conditioning problems 
that may arise in the numerical implementation. A fur- 
ther complication is posed by the general non- convexity 
of the problem, which may result in multiple local min- 
ima that may trap standard descent-based optimization 
techniques (e.g., conjugate gradient 26]) into false so- 
lutions. Therefore, global optimization techniques need 
to be applied, such as genetic [27| . evolutionary [28| . 
or particle-swarm [29| algorithms, whose convergence is 
typically rather slow. Taking into account that each it- 
eration in the optimization procedure may require sev- 
eral numerical solutions of the eigenproblem in (jlip . 
the resulting overall computational burden can become 
prohibitive. From the above considerations, it is clear 
that any blind attempt of attacking such a complex and 
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computationally-demanding problem may be deemed to 
failure. In a step-by-step approach, it appears more rea- 
sonable to start addressing some preliminary issues, such 
as: 

i) A priori estimation of realistic lower bounds for the 
various noise constituents, and comparison with the 
current status and trends, in order to assess the po- 
tential reduction achievable by further optimization 
(and, hence, its worthiness). 

ii) Identification of the effective problem's dimension, 
as a function of the key cavity design parameters. 

Hi) Gathering of prior information (e.g., optimal beam 
profiles and associated structural features) to be 
exploited in order to intelligently drive the opti- 
mization process. 

The rest of the present paper accordingly deals with 
the above issues. 
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FIG. 2: Optimal (minimum-noise) compact-support beam 
profiles in (|21[) . Continuous curve: substrate Brownian (q = 
— 1); dashed curve: coating (q — 0); dotted curve: substrate 
thermoelastic (q = 1). 



III. SOME THEORETICAL BOUNDS 
A. Absolute Bounds: Compact-Support Beams 

1. Variational Problem 

The simplest and crudest approximation of the orig- 
inal optimization problem in (|17[) . leading to a well- 
posed and analytically-treatable problem, consists of as- 
suming the beam profile to exhibit a compact spatial 
support within the [0, 1] interval, thereby implying zero 
diffraction-losses. As we shall see, this implicitly prevents 
the profile from bein| a solution of the eigenproblem in 
(fTTj) . Letting / = \(p\ , one is thus led to the variational 
problem in the space Ll Q ^ of summable functions, 
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{rain) 
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under the constraints 



/: [0,1] 



/ (r) rdr — 1, 



(18) 



(19a) 



(19b) 



whose solution is given below. The arising results are 
anticipated to provide absolute lower bounds, which may 
not be attainable, in view of the mentioned unphysical 
simplifying assumptions. From the Lagrange theory of 
constrained optimization [30| , the constrained variational 
problem in (fT5|) and (|19b[) can be recast into the uncon- 
strained optimization of the Lagrangian functional [40( ] 



A[/,/x]= K*Hi[f\ -2/i 



/ (r) rdr — 1 



, (20) 



where [i is the so-called Lagrange multiplier. It is shown 
in Appendix [X] that this problem admits a unique so- 
lution f s (i.e., an absolute minimum), obtainable using 
variational calculus, viz., 



f s (f) = \^(f)\ 2 = (q + 2)(l-f 2 y. 



1<?<1, (21) 



which also satisfies the positivity constraint in (|19a|) . The 
corresponding (minimum) noise components are given by 



°abs — L 



MI 



l 



(22) 



with r denoting the Gamma (factorial) function [2J, Sec. 
6.1]. The optimal beam profiles are shown in Fig. [21 
whereas the corresponding (minimum) noise values are 
collected in Table [Tl] 

The following remarks are in order: 

• The noise-minimizing beam profiles can exhibit 
step discontinuities, or even singularities at f = 1 
(see Fig. [5]). This is neither surprising (in view 
of the relaxation of the physical feasibility con- 
straints) nor undermining of the mcaningfulncss 
of the preliminary results derived at this stage as 
(anticipated) absolute lower bounds for the actual 
problem. The reader is referred to Sect. IIII Bl below 
for more realistic bounds. 

• For the coating noises (q = 0) , the optimal profile is 
perfectly flat, thereby supporting previous intuitive 
arguments in favor of flat-top beams @, [1] ■ 

• For the substrate noises, the optimal beam profile is 
appreciably rounded (non-flat) for the thermoelas- 
tic component (q = 1). This should be taken into 
account when assessing the performance of config- 
urations featuring sapphire test-masses, for which 
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substrate thermoelastic noise is known to be dom- 
inant llOf - In this framework, use of hyperboloidal 
beams [H], HI, EH as physically-feasible approxi- 
mants should be explored. On the other hand, the 
optimal intensity profile for the Brownian compo- 
nent (q = —1) is close to flat, with a steep increase 
at the mirror's edge. This is clearly unphysical, 
but may be suggestive of using an annular beam. 
Note that the above results pertain to the mini- 
mization of a single noise constituent. Extensions 
to the minimization of a given combination of noise 
constituents are possible, but most likely need to 
be pursued numerically, via suitable discretization 
of the involved operators. 



B. More Realistic (Tighter) Bounds: 
Diffraction-Loss vs. Band-Limitation Constraints 

Besides the diffraction loss constraint, a less obvious 
(and competing) constraint exists, stemming from an ab- 
stract property of the eigenmodes of (III|) : band limita- 
tion. 



1. Band- Limitation Property 

Applying the HT operator at both sides of the eigen- 
problem in (jlip and using the more or less obvious iden- 
tities 

H 1 [f]=n[U(f)m], H[H[f}] = f, (23) 

with II denoting the unit rectangular-window function, 
n(£) = 1, < £ < 1, 11(f) = 0, f > 1, one obtains 

H [<£exp (iV)] (nN D f) = i-^TI (f) exp [-iV (f)] <j> (f) . 

7 

(24) 

Equation (I24[) shows that the HT of the function 
4> exp (tV) (and, a-fortiori, of the function 4>) has a com- 
pact support, vanishing outside [0,ttNd]- Technically, 
the HT plays the role of a wavenumber spectrum, and ac- 
cordingly ttNo is the spatial bandwidth of the field. Note 
that the spatial bandwidth is proportional to the num- 
ber of electromagnetic degrees of freedom Njj in (fT4")l . 
It is therefore natural to try approximating the optimal 
(but, as anticipated, unphysical) beam profiles obtained 
in Sect. IIII A II using a basis in L 2 , with finite spa- 
tial bandwidth ttNd- It is worth stressing that no con- 
structive procedure is given for retrieving a mirror profile 
for which such a superposition is an actual eigenmodc. 
Nonetheless, being a physically admissible (finite spatial 
bandwidth) profile, it is expected to yield tighter noise 
bounds, as compared to (f2"2"]l . 




FIG. 3: PSWF eigenvalues (magnitude) as a function of order 
m, for various values of Nr> ■ The semi-log scale highlights the 
step behavior with exponential tail (see the discussion in Sect. 
IIII B 2|) . Squares: Nd = 1; circles: Nd = 5; up-triangles: 
N D = 10; down-triangles: Nd = 20. 



2. Prolate-Spheroidal Wave-Function Expansion 

A more or less obvious choice for the space-bandlimited 
basis is provided by the so-called prolate- spheroidal wave- 
functions (PSWFs) [U, [32, |33 [H, which satisfy the 
eigenproblem [42| 

fjip(f) = vkNoHx [if] (irN D f). (25) 

In our implementation, the PSWFs are calculated follow- 
ing the approach in [33| . It can be shown that (apart from 
irrelevant complex multiplicative constants) the solutions 
of (|25p . tp n , are real and satisfy the double orthogonality 
condition 

(<Pn,<Pm) = Snmi (<fn, <Pm)i = Vn^nm, (26) 

where S mn is the Kronecker symbol, fj n indicates the n-th 
eigenvalue of (|23|) . and (•, •) and (-, •) 1 denote the L^ Q ^ 
and L^ Q ^ (cylindrical) inner product, respectively. The 
eigenvalue spectrum of (|25p . shown in Fig. [3] for several 
values of Nd, has a step-like behavior: the first ~ Nd 
eigenvalues are close to one in magnitude, while the re- 
maining decay exponentially to zero [32j. The semi-log 
scale utilized in Fig. [3] highlights the step-like behavior 
(with exponentially-decaying tail) of the eigenvalue spec- 
trum. The double-orthogonality condition in ()26|) implies 

/oo 
\^ m {f)\ 2 fdr = \~\n m \ 2 . (27) 

In view of (|2T|) and the noted behavior of the eigenval- 
ues, the first ~ Nd eigenmodes are almost fully localized 
in [0,1], while the remaining ones are almost fully (de- 
vocalized to r > 1. A plot of a few PSWFs of increasing 
order is shown in Fig. [5J for Nd = 11.58. Also shown, 
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FIG. 4: PSWF profiles for Nn = 11.58 and various orders (solid curves), (a): m = 0; (b): m — 2; (c): m = 5; (d): m — 8; (e): 
m — 10; (f): m= 13. Also shown, as reference (dashed curves), is the behavior of the infinite-mirror (GL-type) solutions. 



as a reference, is the behavior of the infinite-mirror (GL- 
type) solutions. It is observed that the agreement be- 
tween the two is rather good for low orders, for which 
the functions are localized, and deteriorates as the order 
m approaches Nd, beyond which the functions exhibit 
the anticipated de-localization. 

The best (in L 2 norm) band-limited approximation 
of the compact-support minimum-noise beam profiles is 
therefore provided by the PSWF expansion 

M T -1 

(t>BL{r) = £ c m^m(r), (28) 
m— 

with the coefficients c m obtained via Fourier-type projec- 
tion m, 



Mt-I 

\ ((^'^«)i) 2 

\ 71=0 

It is readily shown that the truncation order Mt in (|28|) 
is dictated by the prescribed diffraction loss. Under the 
ideal step-like assumption for the eigenvalue dependence 
on index, whereby \fj m \ — 1, Vm < Nd, and fj m = 0, 
Vm > Nd, the diffraction- loss constraint would be sat- 
isfied for any Mt < Nd, however small the prescribed 
Ct- A conservative estimate of Mr, taking into account 
the actual, albeit tiny, departure of the m < Nd eigen- 
value magnitudes from unity may be obtained from the 
obvious inequality 

M T -1 

C^BL] = (1 - \Vm\ 2 )\c m \ 2 

711—0 

M T -1 

< (1- |f7A.^-i| 2 ) J2 \c m \ 2 

m=0 

= (l-fe^-il 2 ), (30) 
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FIG. 5: Truncation index Mt in the PSWF expansion (cf. 
(|32[0 as a function of number of electromagnetic degrees of 
freedom Nn- 

where use has been made of the double-orthogonality 
conditions in (|26|) . the fact that the |?7 m | form a 
monotonically-decreasing sequence, and the unit-norm 
constraint, viz. 

M T -1 

£ \c m \ 2 = HblW 2 = I- (31) 

?TI = 

We accordingly get the following (conservative) estimate 
for the truncation order, which sets the effective dimen- 
sion (number of available design parameters) of the beam 
(mirror) optimization problem: 

M T = largest m : (1 - |?7m-i | 2 ) < £ T - (32) 

For Ct = lppm, the truncation index Mt computed 
from ([50)) is plotted as a function of Nd in Fig. \5\ We 
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FIG. 6: Realistic noise bounds obtained from band-limited 
beam profiles (cf. ([28])) and truncation index Mt as in (|32[) 
(cf. Fig. [5j, as a function of Nd- (a): Substrate Brownian 
(q — —1); (b): Coating (q — 0); (c): Substrate thermoelastic 
(q — 1); Also shown, as references, are the corresponding 
mirror-radius scale (top axis, assuming L = 4km and A = 
1064nm), the absolute bounds (dotted lines, cf. (|22p ). and 
the noise values for a = 16cm (i.e., Nd = 12.03). 



may loosely conclude that the effective dimension of the 
optimization problem is 



(33) 



The inequality in (|30[) will be reasonably tight when rep- 
resenting functions that are essentially localized within 
the unit-disc (mirror-confined beams), whose projection 
onto the de-localized eigenstates with m > Nd is negli- 
gibly small. 

More or less obviously, the accuracy of is strictly 
dependent on (and expected to increase with) the number 
of terms in the truncated expansion. It makes therefore 
sense to check how close one can get to the optimal pro- 
files for various meaningful values of Nd (and hence, via 
(02]), Mr). 

Figures [(Ji-c show the behavior of the noises associ- 
ated with the band- limited profiles in (|28p. referred to as 
Sbl, as a function of Nd- For all three cases, the noise 



FIG. 7: As in Fig. |S] but representative band-limited beam 
profiles for various values of Nd- Continuous curve: Nd = 
6.77 (a = 12cm); dashed: N D = 13.58 (a = 17cm); dotted: 
N D = 22.74 (a = 22cm). 



decreases with increasing Nd , and appears to asymptot- 
ically approach values close to the absolute bounds in 
([2"2"]l . Recalling (fT4"|) , the variation of Nd was obtained 
by tuning the cavity length L and the laser wavelength 
A at the reference values in Adv-LIGO, and taking the 
mirror radius a within a realistic range (see the scale on 
the top axis of Fig. [S]). For each value of a (and, hence, 
Nd), the truncation index Mt (cf. Fig . [5]) was derived 
according to ([32]) , with Ct — lppm |44j |. 

Figure [7] shows the corresponding band-limited beam 
profiles for selected values of Nd- For the q = 
case (coating Brownian and thermoelastic noises) , for in- 
stance, it is observed that, as Nd increases, the profile 
tends to exhibit a more rapid ripple and a steepest decay. 
It can be argued that the flatness of the profile does not 
seem to be an essential ingredient for the coating noise 
reduction. 

To sum up, it is seen that the diffraction-loss constraint 
sets an upper limit to the effective dimension of the opti- 
mization problem, via the finite-spatial bandwidth prop- 
erty of the physically admissible solutions of the cavity 
eigenstate equation. Thus, the only way to approach the 
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absolute minima of the noise constituents acting on the 
beam/mirror profiles, under a prescribed diffraction-loss 
constraint, is by increasing Np, (i.e., if the cavity length 
and laser wavelength are kept fixed, by increasing the 
mirror radius a). However, as seen from Fig. [6j increas- 
ing Njj (aka, a) beyond a certain value pays little, as the 
noise curves roll-off very slowly beyond a certain point, 
and tend to settle. Going, e.g., from a = 16 cm (Adv- 
L1GO baseline design) to a = 23 cm reduces the coating 
noise only by ~ 14%. Moreover, besides the technological 
challenges involved, this raises some model-consistency 
issues related to the actual validity of the underlying 1TM 
approximation. 



3. Validity of the ITM Approximation 

In [2l|, the ITM approximations in (fT]) have been val- 
idated and calibrated against the FTM numerical solu- 
tions in [nj, for GB and MB profiles. Assuming a 40 
kg fused-silica test mass and Cq- — lppm (design speci- 
fications for Adv-LIGO), the ITM predictions for a MB 
profile were found to yield errors < 10% in the coating 
(Brownian as well as thermoelastic) noises and < 25% 
in the substrate Brownian noise (the thermoelastic noise 
component being negligible for fused-silica test-masses 
[Toj|). for mirror radii a < 17cm. For sapphire test- 
masses, the error in the (dominant) substrate thermoe- 
lastic noise component was found to be comparable to 
the substrate Brownian case for fused silica. Assuming 
that comparable figures apply to the band-limited beam 
profiles in (|28p too, some representative values of the re- 
alistic bounds for the case a — 16cm {Nr> = 12.03) are 
reported in Table |n] (scaled to the corresponding abso- 
lute bounds in (|22[)). One observes a moderate increase, 
as compared with the absolute bounds, of a factor ~ 1.14 
for the substrate Brownian noise, ~ 1.31 for the coating 
(Brownian as well as thermoelastic) noises, and ~ 1.55 
for the coating thermoelastic noise. 



IV. COMPARISON WITH CURRENT STATUS 
AND TRENDS: GAUSSIAN AND "MESA" 
BEAMS 



It is suggestive to compare the above derived bounds 
with those attainable by the current (Gaussian) and pro- 
posed (mesa) beam profiles. For these profiles, with- 
out solving the eigenvalue problem in (|11[). one can ex- 
ploit simple approximate analytic solutions for the dom- 
inant eigenmode, valid in the (transversely) infinite mir- 
ror limit, estimating the relevant diffraction losses via the 
so-called "clipping approximation" [34|], i.e., by using the 
infinite-mirror approximate field distributions in the first 
equation in (|16[) . 



A. Gaussian Beams 

The scaled field distribution, in the infinitc-radius- 
mirror approximation, for a GB can be expressed as 



t>GB(r,w ) = Zjqb exp 



2wl 



(34) 



where Eqb is a normalization constant, and the waist 
parameter wq is fixed by the clipping-approximated 
diffraction- loss constraint, 



W c 



(-log£ T ) 



(35) 



In view of the particularly simple analytic expression of 
the field intensity distribution (and of its HT), the scaled 
noise functional in (|10p can be computed in closed form, 



Sgb = 



si n 



\4>. 



GB | 



2*r(§ + i) {-\ogC T f- 
(i-£ T ) 2 



B. Mesa Beams 



> 



-1. (36) 



A MB profile supported by a nearly-flat MH-shaped 
mirror can be synthesized via coherent superposition of 
GBs, with identical waist parameter wq and parallel op- 
tical axes, launched from a circular aperture of radius 
Ro in the waist plane. As shown in [HI, 0], in the 
infinite-mirror approximation, such a beam profile can 
be effectively represented in terms of a Gauss-Laguerre 
(GL) expansion, which, in the scaled form used here, can 
be written as f45| 



4>MB(r, Ro, w , N D ) = E M b exp [iQ (r)] 

OO 

x ^2 A m (R ,wo,N D ) 



m=0 



V2f 



(37) 



In (|37p . Smb is a normalization constant, 0(f) is an 
irrelevant phase distribution, and the expansion coeffi- 
cients A m are given by 



A m {R , w , N D ) = (-l) m P m + 1 



R o 
2wl 



x exp 



2im arctan 



1 



, (38) 



with P denoting an incomplete Gamma function [2J, Eq. 
(6.5.13)]. The orthonormal GL basis functions in (|3~7| 
are 

il> m {£) =V2exp(-^\ L m (e), (39) 
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with L m denoting an mth-order Laguerre polynomial [24l . 
Chap. 22]. In the most general case, the scaled MB pro- 
file depends on three parameters: Rq, wq, and No- How- 
ever, the (clipping-approximated) diffraction-loss con- 
straint introduces a relationship Rq = Rq c (wq , No ) , il- 
lustrated in Fig. [8] for several representative values of 
No, which reduces the number of independent param- 
eters to two (wq,No)- It is worth noticing that, in the 
topical MB literature, the waist parameter u>o is heuristi- 
cally chosen according to a minimum spreading criterion 
EH, viz. 



(MS) 



1 



1 



(40) 



in an attempt of achieving the best tradeoff between top- 
flatness and edge-steepness of the beam intensity profile. 
This further reduces the number of independent parame- 
ters to one (No)- When minimizing the noise functional 
in (|17p . the above choice, while intuitively sound, is not 
necessarily justified a priori from the mathematical view- 
point, and the more general two-parameter optimization 
problem 



=i(mm) 
''MB 



mm 



n 



\4> 



MB 



(41) 



should be considered instead. However, the heuristic 
minimum-spreading criterion in (|40|) is pretty close to 
optimal for sufficiently large No- As an illustrative ex- 
ample, the functional in (|4ip for q — (coating Brownian 
and thermoelastic noises) is shown in Fig. [5] as a function 
of wq (scaled to the reference minimum-spreading value 
in (|4T)|) ), for several representative values of No, within 
the parametric range of potential interest for Adv-LIGO. 
It is observed that the curves at fixed No exhibit a broad 
minimum around (but not exactly at) wq = w MS \ which 
becomes deeper and broader as No is increased. A similar 
behavior is observed for the substrate Brownian (q = —1) 
and thermoelastic (q — 1) noises, and is not shown for 
brevity. 




0.100 0.125 0.150 0.175 0.200 0.225 0.250 0.275 



FIG. 8: Relationship between the MB parameters Ro and wo 
arising from the clipping-approximated diffraction-loss con- 
straint (£t = lppm), for various values of No (assuming 
L = 4km and A = 1064nm). Squares: N D = 6.77 (a = 12cm); 
circles: N D = 9.21 (a = 14cm); up-triangles: N D = 12.03 
(a = 16cm); down-triangles: Nd = 15.23 (a = 18cm); di- 
amonds: N D = 18.80 (a = 20cm); stars: N D = 22.74 
(a = 22cm) . White bullets mark the minimum-spreading con- 
figurations (cf. (TIP)) '). 




C. GB vs. MB vs. Absolute and Realistic Bounds 

The noise levels achievable via a GB profile and a ref- 
erence MB profile (minimum noise, for a = 16cm, i.e., 
No = 12.03, cf. Fig. [9]), scaled to the corresponding ab- 
solute and realistic bounds derived in Sect. IIIH are also 
included in Table |TT] [47] . As already established in [l(| , 
MB profiles yield consistently lower noises than the GB 
counterparts, with reductions of nearly a factor ~ 2.13 
in the coating (Brownian as well as thermoelastic) noise, 
and of ~ 1.45 and ~ 3.07 in the substrate Brownian and 
thermoelastic noise components, respectively. 

By comparison with the absolute and realistic bounds, 
one notes a potential for significant further reductions. 
Specifically, as compared to the MB reference values, the 
realistic bounds indicate potential reductions of nearly 



FIG. 9: As in Fig. [8] but coating (q = 0) noises as a function 
of wo (scaled to its minimum-spreading value) . 

a factor 1.8, 2.5, and 2.9 for the substrate Brownian, 
coating (Brownian as well as thermoelastic) and sub- 
strate thermoelastic noises, respectively, thereby justi- 
fying the further exploration of alternative numerical- 
optimization-driven configurations. 

D. Optimal vs. Good Profiles 

On the basis of the above analysis, a robust (e.g., ge- 
netic) optimization algorithm based on a mirror parame- 
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TABLE II: Comparison between absolute (cf. ([22}) and realistic (for a = 16cm, i.e., No = 12.03, estracted from Fig. [6| noise 
bounds. Also shown, as references, are the noise levels attainable with GB and reference MB (minimum noise, for a = 16cm, 
i.e., N D = 12.03, cf. Fig. [9} profiles. 



q 


b abs 






Q(min) i rf(min) 
MB / S abs 


Sgb/Sbl 


q(min) t q 
°MB 1 °BL 


-i 


1.5708 


1.145 


2.965 


2.043 


2.591 


1.785 





2 


1.313 


6.907 


3.238 


5.256 


2.465 


l 


4.712 


1.552 


13.658 


4.454 


8.801 


2.870 



tcrization consistent with the problem's effective dimen- 
sion, aimed at getting as close as possible to the realistic 
(lower) noise bounds, could be implemented with relative 
ease. 

It should stressed, however, that any optimal design 
should also cope with some more or less obvious addi- 
tional requirements to be also rated as a good design, 
e-g-, 

• The optimal mirror should note pose critical tech- 
nological challenges; 

• The optimal field should be easy to launch, i.e., 
should couple effectively to the injected laser beam; 

• The optimal field should be reasonably stable 
against misalignment and mirror manufacturing 
tolerances. 

Any candidate sub-optimal designs should be ulti- 
mately gauged on the basis of their compliance with the 
above practical requirements. 

V. CONCLUSIONS AND PERSPECTIVES 

In this paper, based on the ITM approximations in 
[20l l2l| , we have addressed some key preliminary issues 
in connection with the optimal beam-shaping problem for 
thermal noise reduction in advanced gravitational wave 
interferometric detectors. The main conclusions can be 
summarized as follows: 

i) The estimated lower-bounds in a realistic config- 
uration, accounting for physical-feasibility-induced 
(diffraction losses, band-limitation) and model- 
consistency (ITM approximation) constraints, in- 
dicate the possibility of significant noise reductions 
(cf. Table |TTJ as compared with the current sta- 
tus and trend. In particular, for the coating noise 
(dominant for the case of fused-silica test-masses), 
a potential reduction of nearly a factor 2.5 is esti- 
mated, as compared with the MB counterpart. 

ii) The key role of the number No of electromagnetic 
degrees of freedom (aka, the Fresnel number) of the 
optical cavity in establishing realistic lower-bounds 
has been highlighted. In this connection, while the 
possibility of increasing No by acting on the cav- 
ity length or the laser wavelength does not appear 



technologically viable for second-generation detec- 
tors (thereby leaving the mirror radius as the only 
adjustable design parameter), it could be taken into 
account for third-generation instruments. 

Hi) From inspection of the band-limited beam profiles 
derived in Sect. IIII B 21 one can infer important 
(sometime counter-intuitive) prior information to 
intelligently drive the optimization process. For in- 
stance, for the coating noises, it clearly emerges 
that the flatness of the beam is not a critical re- 
quirement, since profiles with ripples (cf. Fig. [7]) 
can perform better than flat-top MBs. This obser- 
vation, which is also consistent with the results ob- 
tained using higher-order modes in spherical-mirror 
cavities [13], should be taken into account when 
parameterizing the functional space chosen for the 
optimization problem. 

We believe that the above results pave the way for 
the actual optimization problem, for which a genetic- 
algorithm [27| implementation is currently under inves- 
tigation. Interesting research directions include exten- 
sions of the preliminary study to higher-order (multipo- 
lar, non-axisymmetric) modes. 

While proofreading this paper, we became aware of 
work done independently by M. Bondarescu and Y. B. 
Chen [H, Chap. 3]. Using the GL basis [25[ to synthe- 
size the cavity field, they succeeded in retrieving a spe- 
cial mirror profile which minimizes the coating (Brownian 
and thermoelastic) noise. Remarkably, their results pro- 
vide a nice independent confirmation of the general con- 
clusions drawn here. Indeed, they confirm the key role of 
the mirror radius (aka, cavity Fresnel number) in setting 
the tradeoff between diffraction loss and noise reduction, 
in complete agreement with the general scenario outlined 
here. Also, the minimum noise achieved by their design 
gets pretty close to our corresponding realistic bound. 
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APPENDIX A: MINIMIZATION OF THE 
LAGRANGIAN FUNCTIONAL IN (^0[) 

The stationary solution of the variational problem ob- 
tained equating to zero the functional derivative of (f2"0)) 
can be easily derived using the Gateaux differential 

^ A[/. + eff t /,]-A[/. >M ] =0) WeL i oi] , (A1) 
It is readily verified that 



A[f + eSf,fi] = A[/,/x] + 2e [/] 



— / drrSf (r) 



hi Hi [Sf] . (A2) 



By using (|A2[) in ()A1|) . and interchanging the spectral 
and spatial (HT) integrals in the inner product, the sta- 
tionarity condition can be equivalently written as 



dr 



dRR q+1 J (Rf) 



x / df'f'f s (f')J (Rf')- fi 



rSf(f) = 0, 



W e i[o,i]. 



(A3) 



from which it follows that the stationary profile f s satis- 
fies the integral equation 



dRR q+1 J (Rf) I dr'f' f s (f')J (Rf') = fi, < r < 1. 
i) Jo 

(A4) 

Equation (|A4[) can be solved in closed-form, by inspec- 
tion. Indeed, we capitalize on its nested-HT structure to 
use the following integral identities 

I K*J 5+1 (K)Jo(fK)dft = 2*r(| + l), 

< f < 1, - 1 < 5 < 1, (A5) 



(see, e.g., [13, Eq. (2.12.31.1)]), and 



K 2 J| + l(«): 



dfV(l -f /2 )ij (Rf'), 



q > -1, (A6) 



(ibid, [33, Eq. (2.12.3.6)]). By combining JSgJ) and ([M 
one obtains 
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[r- 



f + i) 



dRR q+1 J Q (Rf) 



x / rff'f'(l - f' 2 )^ Jo(Kf') = 1, 
Jo 

< f < 1, 



which, by comparison with (j A4[) . yields 

(1 -f 2 )" 



AW = m 



2i [T(§ + 1] 



, < f < 1. 



(A7) 



(A8) 



The as yet unspecified multiplier in (|A8[) is chosen so 
as to satisfy the normalization constraint in (|19b[) . 



r ( | + i 



thereby yielding the final result in f^J 
From (TA~2l) and dA4l). it follows that 



(A9) 



S[g]-S[f s }= R*H x \g-f s ] 



V.9 S Lf 0A] : / g(r)rdr = 1. 



(A10) 



Therefore, from the readily verifiable positive- 
definiteness of the functional norm in (|A10j) . one 
concludes that that the stationary profile f s in (|2"Tj) 
yields the absolute minimum of the noise functional in 

na. 
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This investigation does not deal with tilt (in)stability is- 
sues, for which nearly-flat and nearly-concentric configu- 
rations exhibit significantly different responses 
The concept of electromagnetic degrees of freedom [z| is 
related to the minimum number of basis functions that 
are needed to represent within a bounded region, with 
a prescribed error in a given metric, an electromagnetic 
field produced by finite-support (impressed or equivalent) 



The satisfaction of the positivity constraint in (|19a[) . not 
explicitly enforced, is verified a posteriori. 
Throughout the paper, this more generic denomination 
will be utilized, instead of the more detailed "PSWFs on 
a disc" in [H. 

The PSWF basis is the spectrum of the eigenproblem in 
(|11|) with V = 0, representing the eigenmodes of a confo- 
cal spherical Fabry-Perot cavity with finite-size mirrors. 
In our approach, based on the PSWF computation via 
the polynomial expansions in [33|. the scalar products 
in (|29|) can be efficiently computed in a semi-analytic 
fashion. 

In view of the inherent discreteness of the problem, it 
was not always possible to find values of Mt so that 
the diffraction-loss constraint was tightly satisfied (i.e., 
diffraction loss smaller but close to the targeted value). 
In order to ensure the comparison among configurations 
yielding comparable diffraction losses, we slightly relaxed 
the constraint, looking for values of Mt yielding diffrac- 
tion losses closest (~ ±0.3ppm) to the targeted value. 
A similar GL expansion was also derived for analytic 
parameterization of the MH-shaped mirror-deformation 
profile (cf. [H, Eq. (34)]). 

The Rayleigh distance of the GBs launched is chosen as 
coincident with the distance from the aperture plane to 
the mirror. 

Note that, although the GB and MB results in Table 
iHl have been obtained, for computational convenience, 
using the infinite-mirror approximation in (I37|l , they were 
checked to agree to within 15% with those obtained from 
the numerical solution (via the Nystrom method [35|]) of 
the eigenproblem in pip , for a number of representative 



